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Abstract —We consider a cooperative relaying system which 
consists of a number of source terminals, one shared relay, 
and a common destination with multi-packet reception (MPR) 
capability. In this paper, we study the stability and delay 
analysis for two cooperative relaying schemes; the sensing-based 
cooperative (SBC) scheme and the decision-based cooperative 
(DBC) scheme. In the SBC scheme, the relay senses the channel 
at the beginning of each time slot. In the idle time slots, the 
relay transmits the packet at the head of its queue, while in the 
busy one, the relay decides either to transmit simultaneously with 
the source terminal or to listen to the source transmission. The 
SBC scheme is a novel paradigm that utilizes the spectrum more 
efficiently than the other cooperative schemes because the relay 
not only exploits the idle time slots, but also has the capability 
to mildly interfere with the source terminal. On the other hand, 
in the DBC scheme, the relay does not sense the channel and 
it decides either to transmit or to listen according to certain 
probabilities. Numerical results reveal that the two proposed 
schemes outperform existing cooperative schemes that restrict 
the relay to send only in the idle time slots. Moreover, we show 
how the MPR capability at the destination can compensate for the 
sensing need at the relay, i.e., the DBC scheme achieves almost 
the same stability region as that of the SBC scheme. Furthermore, 
we derive the condition under which the two proposed schemes 
achieve the same maximum stable throughput. 

I. Introduction 

U SER demands for high data rates and services are 
expected to increase exponentially in the next decade. 
According to the Federal Communication Commission (FCC), 
about 70% of the allocated spectrum in the US is not efficiently 
utilized. Hence, effective utilization of the available spectrum 
is a critical issue that has recently gained great attention. 
Cognitive radios and cooperative diversity have emerged as 
promising techniques to improve the wireless network perfor¬ 
mance in an attempt to exploit the unutilized spectrum ifT)- 
0. Intuitively, by relaying the messages and emptying the 
queues of the primary sources, the secondary node creates 
more opportunities for its own transmission. 

Cooperative diversity is a new paradigm for wireless net¬ 
works, and hence, deep investigation is needed to fully under¬ 
stand the impact of this new paradigm on different network 
layers. Most of the work on cooperative communication has 
focused on the physical layer aspects of the problem 01 • Other 


works El- 0, however, have implemented cooperation at 
the network protocol level, and performance gains in terms 
of stable throughput, average delay, and energy efficiency 
were illustrated. In EL the authors proposed a novel cognitive 
multiple-access strategy, where the relay exploits the bursty 
nature of the transmission of the source terminals via utilizing 
their periods of silence to enable cooperation. In this strategy, 
no extra channel resources are allocated for cooperation, 
and hence, this improve the spectral efficiency. Although 
the proposed strategy provides significant performance gain 
over conventional relaying strategies, the relay is restricted 
to transmit only in the idle time slots. However, allowing 
the relay to send simultaneously with the source with certain 
probability can improve the network performance. 

Utilization of multi-packet reception (MPR) capability has 
received considerable attention in the literature. An MPR 
model was first introduced in 0. Multi-access channel (MAC) 
systems with MPR capability have been addressed in the 
literature in different contexts m- a. However, most of 
which do not deal with cognitive or cooperative systems. In 
0, a MAC network with two primary users, a cognitive relay, 
and a common destination was considered with a symmetric 
configuration. The primary users, simultaneously, access the 
channel to deliver their packets to a common destination, i.e., 
the relay and the destination have MPR capability. The authors 
assume that the relay perfectly senses the channel, i.e., it 
transmits during idle time slots, where the primary nodes are 
not transmitting. In this scheme, the authors assume that the 
primary users transmit simultaneously in each time slot which 
may cause degradation in the network performance especially 
in the presence of weak channels. Moreover, as the number 
of primary users increases, the complexity of the relay and 
the destination increases because these nodes have to decode 
the message of all the transmitting nodes in each time slots. 
In fTOll . the performance of an Ad-Hoc secondary network 
with N secondary nodes accessing the spectrum licensed to a 
primary node was demonstrated. Both cases of perfect and 
imperfect sensing were considered. In the perfect sensing 
case, the secondary nodes do not interfere with the primary 
node and thus do not affect its stable throughput. However, 
with imperfect sensing, the secondary nodes control their 
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transmission parameters, such as the power and the channel 
access probabilities, to limit the interference on the primary 
node. To compensate for the effect of interference, the authors 
explore the use of the secondary nodes as relays for the 
primary node traffic. 

The average delay encountered by the packets is one of 
the most important metrics in evaluating the performance of 
wireless networks. In hd, the delay analysis for a network 
consists of one primary user, one cognitive relay, and a com¬ 
mon destination was presented, where the relay is restricted 
to send only in the idle time slots with full priority for 
the relaying queue. Moreover, in fl2l . the delay analysis for 
the same network with a randomized cooperative policy was 
investigated, where the relay node serves either its own data 
or the primary packets with certain service probabilities. The 
authors, in tm, showed that the randomized policy enhances 
the cognitive relay delay at the expense of a slight degradation 
in the primary user one. 

In this paper, we consider a general number of source 
terminals, one half duplex relay, and a common destination 
that has MPR capability, i.e., the destination can decode the 
message of more than one transmitting node in the same 
time slot. We consider a slotted time division multiple access 
(TDMA) framework in which each time slot is assigned 
to one source terminal only. We propose two cooperative 
schemes. The first scheme is the sensing-based cooperative 
(SBC) scheme, where the relay senses the channel at the 
beginning of each time slot. If the relay detects an idle 
time slot, the relay transmits the packet at the head of its 
queue. Alternatively, if the relay detects a busy time slot, it 
decides probabilistically to either listen to the source packet 
and store it if the destination fails to decode it successfully or 
to interfere with the source terminal transmission. We optimize 
this probabilistic scheme to maximize the network aggregate 
throughput and characterize the stability region. The main 
difference between this scheme and that in S3 is that, in El, 
the authors restricted the relay to send only in the idle time 
slot. Moreover, unlike 0, the relay interferes with the source 
terminals with certain probability to limit the adverse effects 
of the interference. 

In the SBC scheme, the relay depends on the sensing 
information to decide either to listen or to transmit. Although 
we do not take into consideration the sensing errors and the 
consumption of power, in practice these factors can negatively 
affect the performance of the system. Hence, we propose 
the decision-based cooperative (DBC) scheme, which is the 
second cooperative scheme proposed in this work. In the 
DBC scheme, unlike ED that assumes imperfect sensing, the 
relay does not sense the channel and it decides, according to 
certain probabilities, either to listen to the source or to transmit 
whether the time slot is idle or busy. It is worth noting that 
the complexity of our proposed schemes does not increase as 
the number of source terminals increases, unlike 0, because 
for any number of source terminals the destination decodes, at 
most, the packets of two transmitting nodes; one of the source 
terminals and the relay. 

In this work, we focus on the medium-access layer and 
address the impact of the proposed schemes on multiple-access 


performance metrics such as the stable throughput region and 
average delay. We show that in the two proposed schemes, 
the queues of the source terminals and those of the relay 
are interacting. Since the stability analysis for more than 
two interacting queues is difficult, we resort to a stochastic 
dominance approach. The stability analysis of interacting 
queues was initially addressed in |[T3l . and later in f]~4l . where 
the dominant system approach was explicitly introduced. The 
average delay is also an important performance measure, 
and its analysis illustrates the fundamental trade-off between 
the rate and reliability of communication. Delay analysis for 
interacting queues is a notoriously hard problem that has 
been investigated in lfl5l and in 0 for ALOHA with MPR 
channels. Hence, we also utilize stochastic dominance to 
approximate the average delay of the proposed schemes. 

Our contributions in this paper can be summarized in the 
following points 

• For the two proposed schemes, we present the stability 
analysis and derive the stability conditions for each queue 
in the system. We formulate an optimization problem 
to maximize the weighted aggregate stable throughput 
of the network and characterize the stability region via 
optimizing the probability of each action taken by the 
relay. The problem is formulated as non-convex quadratic 
constrained quadratic programming (QCQP) optimiza¬ 
tion problem m. We use the feasible point pursuit- 
successive convex approximation (FPP-SCA) algorithm 
EO to achieve a good feasible solution by approximating 
the non-convex constraints as linear ones. 

« We analyze the average delay performance for the two 
proposed schemes and derive approximate delay expres¬ 
sions using the dominant system approach. 

• We show that the SBC scheme provides significantly 
better performance over existing cooperative schemes as 
in 0- Moreover, the SBC scheme exploits the unutilized 
spectrum more efficiently than other schemes, because 
the relay not only transmits in the idle time slots but 
also has the capability to, simultaneously, transmit its 
packets with the source terminals. The relay uses this new 
attribute in a mild way to mitigate the negative effects of 
the interference and to enhance the maximum aggregate 
stable throughput of the network. 

• We demonstrate that the DBC scheme, in certain cases, 
achieves the same stability region achieved by the SBC 
scheme. We also illustrate how the MPR capability at 
the destination can compensate for the need for the 
relay to detect the idle time slots to transmit its packets. 
Furthermore, we show that removing the MPR capability 
from the destination, in absence of sensing at the relay, 
causes catastrophic degradation in the performance of the 
system. 

« We derive the channel condition under which the two 
proposed schemes achieve the same maximum stable 
throughput. Under this condition, sensing the channel by 
the relay becomes useless, and hence, the two proposed 
schemes provide the same performance. 

The remainder of the paper is organized as follows. In 
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Section El we describe the system model. The SBC scheme 
with its stability and delay analysis is introduced in Section 
Hill followed by the DBC scheme in Section UV] Numerical 
results are then presented in Section [Vj We demonstrate how 
the MPR capability at the destination can compensate for the 
relay need to sense the channel in Section [VI] Finally, the 
paper is concluded in Section 1 VIII 

II. System Model 

We consider the uplink of a TDMA system that consists of 
M source terminals {si}^£ 1 , one shared relay (r), and one 
common destination ( d ), as shown in Fig. [T] We define the set 
of the transmitting nodes T={S, r}, where S={si, sm} 
is the set of source terminals, and that of the receiving nodes 
L={r, d). The source terminals access the channel by dividing 
the available resources among them, i.e., time slots in this 
case. Each terminal is allocated a fraction of the time. Let w.i 
denote the fraction of time allocated to the source terminal s t , 

where i £ {1,2,., M}. We assume continuous values of the 

resource sharing vector w = [wi,W 2 , Hence, we can 

define the set of all feasible resource sharing vector as follows 

M 

A = jut = [w 1 ,W2,...,w M ] G R+ , Wj = 1 j. (1) 

i—1 

First, we describe the physical layer model. All wireless 
links are assumed to be stationary, frequency non-selective, 
and Rayleigh block fading. The fading coefficients h m , n , 
where m£T and n£L, are assumed to be constant during 
each slot duration, but change independently from one time 
slot to another according to a circularly symmetric complex 
Gaussian distribution with zero mean and variance p 1 tn n . All 
wireless links are corrupted by additive white Gaussian noise 
(AWGN) with zero mean and unit variance. All nodes transmit 
with fixed power P. An outage occurs when the instantaneous 
capacity of the link (to, n) is lower than the transmission rate 
R. Each link is characterized by the probability 

fmn = P{I? < log 2 (l+P|/l mjn | 2 )} = exp ( ~^p~2 - ) ( 2 ) 

\ * Pm,n / 


which denotes the probability that the link (to, n) is not in 
outage. Let g denote the probability that the link (to, n) is 
not in outage in presence of interference from node /, where 
to, I € T and n £ L. In 0, the authors derived the same term 
but for symmetric configuration. We relax this assumption and 
re-derive the term to fit our model in Appendix A. 

Next, we describe the medium access control layer model. 
Time is slotted with fixed slot duration, and the transmission of 
a packet takes exactly one time slot. Each source terminal has 
an infinite buffer (queue) to store its own incoming packets. 
Packet arrivals of the source terminals are independent and 
stationary Bernoulli processes with means A, (packets per 
slot), where i £ {1,2,....,M}. The relay has M relaying 
queues (Q ri , Qr 2 , • • ■> Qr M ) to store the packets of the source 
terminals that are not successfully decoded at the destination. 
Let Q\ denote the number of packets in the l-th queue at the 
beginning of time slot t. The instantaneous evolution of the 
l-th queue length is given by 

Q* +1 = (Qi-Yi) + + xt ( 3 ) 

where l £ {si, r \,..., tm} and (a:) + = max{.T, 0}. The 
binary random variables Yf and Xf, denote the departures and 
arrivals of Qi in time slot t, respectively, and their values are 
either 0 or 1. 

III. Sensing-Based Cooperative Scheme 

In this section, we introduce the proposed SBC scheme 
followed by its stability and delay analysis. We assume that 
the relay can sense the channel, and perfectly detect the idle 
time slots. Moreover, we assume that the errors and delay 
in packet acknowledgement (ACK) are negligible, which is 
reasonable for short length ACK packets as low rate codes can 
be employed in the feedback channel (4). In the SBC scheme, 
the source terminals operate according to the following rules: 

• Each terminal transmits the packet at the head of its queue 
in its assigned time slot, whenever the queue is not empty. 

« If the destination decodes the packet successfully, it sends 
an ACK which can be heard by both the transmitting 
terminal and the relay. The terminal drops this packet 
upon hearing the ACK. 

• If the destination does not receive the packet successfully 
but the relay does, then the relay stores this packet at the 
corresponding relaying queue and sends an ACK to the 
source terminal. Afterwards, the relay is responsible for 
conveying this packet to the destination. 

• If a packet is received successfully by either the destina¬ 
tion or the relay, the packet is removed from the termi¬ 
nal’s queue. Otherwise, the source terminal retransmits 
the packet in its next assigned time slot. 

It is worth noting that if the relay and the destination decode 
a packet successfully, the relay does not store this packet 
in its queue because the packet is already delivered to the 
destination. 

Next, we illustrate the transmission and reception policy of 
the relay using the SBC scheme. Let us consider the time slots 
allocated to the source terminal Si. The relay takes one of the 
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following two actions after sensing the channel to detect the 
queue state of Q Si . 

• If Q Si is empty, i.e. s, is not transmitting, the relay 
transmits the packet at the head of Q ri . 

• If Q Si is not empty, the relay transmits a packet from Q rj , 
where j £ {1,2,M}, with probability 8i] or listens to 

M 

source terminal transmission with probability 1— E 8 tJ . 

3 =i 

Hereafter are some important remarks. First, the relay stores 
the packet of the source terminal if it decides to listen to the 
terminal’s transmission and successfully decodes the packet 
while the destination fails to decode it. Second, we assume 
that the relay is half duplex, and hence, it can not transmit 
and receive a packet in the same time slot. Thus, the relay 
can not receive packets when it decides to interfere with the 
source terminal. 

From the given description of the proposed cooperative 
scheme, it is clear that the decision taken by the relay depends 
on the queue state of the source terminals, and this causes 
an interaction between the relay’s queues and those of the 
source terminals. Stability of interacting queues is complex 
E), thus, we resort to using the dominant system approach. 
To perform the stability analysis of s- tl we assume that all 
other source terminals Sj, where j i, are saturated, i.e, 
their queues always have packets. Moreover, we assume that 
Q ri transmits a dummy packet whenever the relay decides 
to interfere with the source transmission and Q ri is empty, 
where i £ {1,2This dominant system simplifies the 
stability analysis and provides an inner bound on the stability 
region. 


A. Stability Region Analysis 

In this part, we characterize the stability region for the 
SBC scheme taking in consideration the dominant system 
described above. Then, we formulate an optimization problem 
to maximize the aggregate throughput and characterize the 
stability region via optimizing the values of {/3subject 
to constraints that ensure the stability of the system. The 
stability of the overall system requires the stability of each 
individual queue. From the definition in m, the queue is 
stable if 


lim P {Q* < x} = F(x) and lim F(x) = 1 (4) 

t —^oo x—»oo 

We can apply Loynes’ theorem to check the stability of a 
queue HD- Loynes’ theorem states that if the arrival process 
and the service process of a queue are strictly stationary, then 
the queue is stable if and only if the average service rate is 
greater than the average arrival rate of the queue. 

In the dominant system, a packet departs Q Si , where i £ 
{1,2,...,M} in two cases. First, if it is successfully decoded 
by the destination when the relay decides to interfere with Sj 
from one of its relaying queues. Second, if it is successfully 
decoded by at least one node, i.e., the destination or the relay 
when the relay decides to listen to s,. Thus, the average service 


rate of Q Si for certain w L is given by 


M M 

l^i= w i(Yl Pij9 r s id + (! - ftjX/sid + (! - fs,d)f Si r)J ■ 
3 =1 j -1 

(5) 

For stability of Q Si , the following condition must be satisfied 

M M 

A ^9 r Si d + (! - Pij)(f»id + (! “ fsid)fsir)^ ■ 

3 =1 3 =1 

( 6 ) 

A packet arrives at Q ri if the following two conditions are 
met. First, an outage occurs in the link between s,; and the 
destination node while no outage occurs in the link between s, 
and the relay and this happens with probability (1 —fs id )fsir- 

Second, Q Si is not empty and the relay decides to receive from 

M 

Si and this happens with probability — (l — ^ 8 ^). Thus, the 

^ j =l 

average arrival rate of Q ri is given by 

A M 

X ri = W l (l-f Sid )fs i r—{l ~ y^Pij) (7) 

Mi 3=1 

In the dominant system, the packet departs Q ri in three 
cases. First, when Q Si is empty a packet departs Q r , if no 
outage occurs in the link between the relay and the destination. 
Second, when Q Si is not empty a packet departs Q ri if the 
relay decides to interfere with s t with a packet from Qr, - 
Third, if the relay decides to interfere with saturated source 
terminal Sj, where j ^ i, with a packet from Q ri . Thus, the 
average service rate of Q ri can be expressed by 

H ri =Wi((l-— )frd+—Pii9r^]+ W 3p3i9rd- ( 8 ) 

V IM J 3^,3=1 

For the stability of Q ri , the service rate must be higher than 
the arrival rate, i.e., X ri < /i ri , and hence, we have 


Wifrd + E WjPjidrd 
3=1,jAi 
M 

1 - E Aj)(! - f Si d)f Si r + frd - Pii9 S r d ) 

3=1 J 

(9) 

Let A* denote the maximum stable throughput for the i- 
th source terminal at a certain value of Wi. To guarantee the 
stability of the entire network the conditions in ([6ji and ([9} 
must be satisfied. Hence, the maximum stable throughput is 
given by 

A* < min{m,n Ui } (10) 



where can be obtained easily from © as follows 


— 


M 

mfrd + E "ll-ij'lir’d 

3=1, iAi 


, M \ 

w i ((! - E Pij ){ 1 - fSid) fSir + frd ~ Pii 9 rd) 


(ID 


Next, we optimize the values of {/E'}fj=i to achieve 
the maximum weighted aggregate throughput of the network 
subject to constraints that ensure the stability of all queues. 
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Thus, the optimization problem can writterQ. for certain wt, 
as follows 


max 

{/3>5=i 

subject to 


M 

2=1 

M 

^ X > 

3 = 1 

A s 5: Mil 
A s ^ M«s > 

0 < fa < 1, 


e {l, 2,M} 


where a;, is the weight assigned to A*. To obtain the sta¬ 
bility region, we have to vary the value of each x t , where 
i £ {1, 2,..., M}, from zero to one for each w-i. Then vary Wi 
also from zero to one to scan the whole stability region, i.e., 
the convex hull of all the obtained values of {A* }^L 1 provides 
the stability region. On the other hand, we set all values of x. t 
by ones to calculate the maximum aggregate stable throughput 
of the network. 

In the above problem, the summation in the first con¬ 
straint is the probability that the relay transmits a packet 
from one of its relaying queue while Q Si is not empty. The 
second and the third constraints guarantee the stability of the 
queues in the network. In order to solve this optimization 
problem, we define a new M(M + 1)-dimensional vector 
ft=[ftii ,..., /3im, •••) 0mi, /3mm, Aj, A^] T and rewrite 
the optimization problem, in the standard form |20l , as follows 


min x T 3 
P 

s.t. 0 < bj (3 < 1, i £ 

vfp + Ui<0, i £ {1, 2,..., M} (!3) 

(3 T Ai[3 + C J(3 + di < 0, i £ {1, 2,..., M} 

0 A (3 < 1 


where x=[0, 0,..., 0, — xi, ..., —Xm], and the terms i >j, m, Ai, 
Ci, and di can be easily obtained from © and ifTTb . From 
©, the stability of Q Si is represented by linear constraints 
in (3 while, from CD, the stability of Q ri is represented 
by quadratic constraints in (3. The objective and the linear 
constraints are convex. However, the quadratic constraints are 
not convex because A, is an indefinite matrix. In general, 
non-convex QCQP problems are NP hard ED, except for 
special cases such as those in El . Several methods have 
been proposed to approximate non-convex QCQP problems, 
including semi-definite relaxation (SDR) H21 1 . the reformula¬ 
tion linearization technique (RLT) ]23l , and successive convex 
approximation (SCA) l24l . In our case, we use an iterative 
algorithm to obtain a good feasible solution as in Ifl7l . We 
approximate the feasible region through a linear restriction of 
the non-convex parts of the constraints. The solution of the 
resulting optimization problem is then used to compute a new 
linearization and the procedure is repeated until convergence. 
Using the eigenvalue decomposition, the matrix A, can be 


expressed as A ?; = Af + A i , where Af y 0 and A i < 0. 
For any y £ ^M(M+i)xi, we can replace the non-convex 
constraint in (fl3l > by the following convex one 

P TA t(3 + 2y T A~[3 + cf 13 + di < y T A~y (14) 

See El for more details. Thus, the non-convex problem is 
converted to a convex one, and we use Algorithm 1 to solve 
the optimization problem in (El- 


Algorithm 1: We use the FPP-SCA algorithm to achieve 
good feasible point 

t Initialization: set k = 0 and y 0 = 0. 

2 Repeat 

1) solve 

min x T (3 
P 

s.t. 0 < bf (3 < 1 
vj(3 + m < o 

/3 A t ft + 2 yl Aj (3 + cj(3 + di < yj. A i y k 
0 < ft < 1, i£ {1,2 

2) Let [3* k denote the optimal (3 obtained at the k-th 
iteration, and set y k+1 = (3* k 

3) Set k = k + 1. 

until convergence. 


It is worth noting that we can use other techniques to solve 
this problem, but the FFP-SCA obtains good feasible point 
even for very large M. Moreover, a few iterations are required 
for convergence 03. 


11. Average Delay Characterization 

In this subsection, we characterize the average delay for 
the dominant system of the SBC scheme. We derive an 
approximate expression of the delay as in 0, where we 
assume that the relay queues are discrete-time M/M/1 queues. 
If the packet is directly transmitted from the source terminal 
to the destination, it experiences a queueing delay only in the 
terminal’s queue. Alternatively, if the packet is delivered to 
the destination through the relay, this packet experiences two 
queuing delay; one in the terminal’s queue and the other in 
the relay’s queue. The packet experiences only the queueing 
delay at the source terminal with the following probability 

, M s M 

(l~ S Pijjfsid + Yl ftijg r Si d 

= Wi - 3 — - 3 — - (15) 

Mi 

which is the probability that the packet is successfully decoded 
by the destination given that it is dropped from the source 
terminal. Thus, the average delay encountered by the packets 
of the i-th source terminal is given by 


1 We do not restrict the second and the third constraints to be satisfied with 
strict inequality as in {6} and HJ. However, the solver uses the interior point 
method which provide a strictly feasible point. 


Di — eiT Si + (1 — ei)(T Si + T ri ) 

= T Si + (1 - a)T ri (16) 
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where T Si and T ri denote the average queueing delays at Sj 
and r-i, respectively. Since the arrival rates at Q Si and Q ri are 
given by A * and \ ri , respectively, then applying Little’s law 
yields 

T Si = Ni/X h T ri =N r JX ri (17) 

where N. t and N ri denote the average queue size of Q Si 
and Q ri , respectively. The queues of the source terminals 
are discrete-time M/M/1 queues with Bernoulli arrivals and 
Geometrically distributed service rates, and we assume that 
Q ri is a discrete-time M/M/1 queue. Thus, we can easily 
calculate N Si and N ri , by applying the Pollaczek-Khinchine 
formula 1 125 1 , as follows 


Ni = 


~X? + A* 


N r , = 


—At! + A r 


( 18 ) 


Hi X i Hr, A r 

Substituting (ITTb and (IT8] > in (fl6l >. we can write the average 
queueing delay for the i-th source terminal as 

1 A/ .. ,1 Ar- 


IX, = 


Hi X i 


+ (1 - Ci) 


jj.x r 


(19) 


where the expressions of Hi and X ri , and Hn are obtained in 
(0, 0, and ®, respectively. 


IV. Decision-Based Cooperative Scheme 

In this section, we present the proposed DBC scheme 
together with its stability and delay analysis. In this scheme, 
the behavior of the source terminals is exactly the same as in 
the SBC scheme, where each source terminal transmits only in 
its assigned time slots. The source drops the packet if it hears 
an ACK from the relay or the destination, otherwise, the source 
retransmits the packet in the next assigned time slot. The main 
difference between the DBC and the SBC schemes is in the 
way the relay operates. In the SBC scheme, the relay decides 
its operation policy depending on the sensing information. 
However, in the DBC scheme, the relay does not sense the 
channel and it decides randomly either to transmit or to listen 
regardless the queue state of s,. To illustrate the operation 
policy of the relay in the DBC scheme, let us consider the 
time slots allocated to .s, where the relay takes one of the 
following actions: 

• The relay transmits a packet from Q rj with probability 
atij, where j G {1,2 

• The relay listens to the source transmission with prob- 

M 

ability 1 — X] The relay stores the source’s packet 
3—i 

if the relay successfully decodes it while the destination 
fails. 

When the relay decides to transmit without listening to the 
channel, the probability of successful transmission depends 
on the queue state of the source terminal because if the 
terminal’s queue is empty, there is a higher probability that 
the destination decodes the relay’s packet. Hence, it is clear 
that there is an interaction between the relay queues and 
those of the source terminals. To simplify the analysis, we 
assume the same dominant system used in the SBC scheme, 
where Q ri transmits a dummy packet whenever the relay, 
randomly, decides to transmit a packet from Q ri while this 


queue is empty. Furthermore, to perform the stability analysis 
of a certain source terminal, we assume that all other source 
terminals are saturated, i.e., their queues are not empty at 
any time slot. This dominant system simplifies the stability 
analysis and provides an inner bound on the stability region. 


A. Stability Region Analysis 

In this subsection, we follow the same steps as in the 
SBC scheme. First, we characterize the stability conditions 
for all the queues in the DBC scheme. Then, we formulate 
an optimization problem to maximize the weighted aggregate 
throughput and characterize the stability region by optimizing 
the values of under the stability constraints. 

In the dominant system, a packet departs Q Si , where i G 

{1, 2,., M}, in two cases. First, if it is successfully decoded 

by at least one node, i.e., the destination or the relay when 
the relay decides to listen to s,;. Second, if it is successfully 
decoded by the destination when the relay decides to interfere 
with Sj. Thus, the average service rate of Q Si for certain w, t 
is given by 


M M 

Hi = Wi ^(1 — "y ' Q-ij ) (fSid + (1 ~~ fSid) fSir) + ^ ' ot -ij9s i d'j 
3 =1 3=1 

( 20 ) 

For stability of Q Si , the following condition must be satisfied 

M M 

X i < tUi ^(1 — y ' a ij)(fsid + (1 — fSid) fSir) + ^ ^ a ij9s i d'j 

3=1 3=1 

( 21 ) 

A packet arrives at Q ri if the following conditions are 
met. First, an outage occurs in the link between s,; and the 
destination node while no outage occurs in the link between Si 
and the relay. Second, Q Si is not empty which has a probability 
of Xi/Hi, and the relay decides to listen to Si which happens 

M 

with probability 1 — X] a ij- Thus, the average arrival rate of 

3=1 

Q ri is given by 

A M 

X ri = Wi(l - f Si d)f Si r — (1 - a x) (22) 

3=1 

In the dominant system, a packet departs from Q ri in two 
cases. First, if the relay randomly decides to transmits a packet 
from Q ri on the time slot allocated to Si. In this case the packet 
departs Q ri with probability Wi(auf r d{l ~ ^)+a ii g s r ) 1 ^)■ 
Second, if the relay decides to transmits a packet from Q ri 
in the time slot allocated to the saturated source terminal 
Sj, where j X i. In this case the packet departs Q ri with 

M 

probability X! w j a 3i9rd- Thus, the average service rate 

3 = 1 , 3 #* 

of Q ri can expressed by 



\ | Si 

-)+a»£Ld- 

An An 


M 

X! W 3 a 3i9rd 
3=1,3#* 

( 23 ) 
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For the stability of Q ri , the service rate must be higher than 
the arrival rate, i.e., X ri < fi ri , and hence, we have 

m 

WiOiiifrd+ Y w i a ji9rd 

Aj < -—- Ui 

, M S r 

d'j ( (1 Yj ^ij) fSir( 1 fSid) T &iifrd ^^9 r dj 

V j=l 7 

(24) 

To guarantee the stability of the network the following 
condition must be satisfied 

A* < nun{/jbi,/jb Ui } (25) 

where \i Ui obtained directly from (f24l i as follows 

m 

V Ui = Z M ~ fii 

tUjf (1 Yj ^ij) fSir{ 1 fSid) Qiifrd ^iidrd) 

(26) 

We can write the optimization problem to calculate the 
maximum weighted stable throughput of the network as 

M 

max > Xi\\ 

M 

subject to 0 < E a ij * £ {1> 2).> M} 

3 =1 

a; < l-H, is {1,2, . ,M} 

K < Muo * e {1,2,. ,M} 

0 < a Kj < 1, *,j€{l,2, ., M} 

(27) 

This optimization problem is exactly the same problem as that 
in the SBC scheme. Thus, we follow the same steps to solve 
this problem, and use the FPP-SCA algorithm to obtain good 
feasible solution. 


B. Average Delay Characterization 

In the DBC scheme, the source packets are either transmit¬ 
ted directly to the destination or through the relay. We define 
T,, which is the probability that the packet is successfully 
decoded by the destination given that it is dropped from the 
Si, as 

, M N M 

(l~ Y a ij)fsid + Y a ij9 r Si d 

Ti = Wi -—-—- (28) 

Hi 

Thus, the packets of the i-th source terminal experience the 
following average delay 

Di = T Si + (1 - Ti)T ri (29) 

As in the SBC scheme, we can calculate T Si and T ri by as¬ 
suming that Q ri is a discrete-time M/M/1 queue and applying 
Pollaczek-Khinchine formula. Hence, we can write the average 
queueing delay for the i-th source terminal’s packets as 

A = —~ t : + (i - r i)——-w- (30) 

fli Ai 

where the expressions of /x*, \ ri , and /.x rj are obtained in < 12()b . 
(l22l i. and (f23l) . respectively. 


V. Numerical Results 

In this section, we evaluate the performance of the two pro¬ 
posed schemes for the case of two source terminals, i.e., M= 2. 
We compare the performance of the SBC and DBC schemes 
with other cooperative schemes. First, we compare with the 
cooperative cognitive multiple-access (CCMA) scheme, pro¬ 
posed in 11, where the authors assume that the relay only 
transmits in the idle time slots. Moreover, the relay only helps 
the terminals which, on average, have worse channel condition 
than the relay itself. In other words, the relay assists the 
terminals whose outage probability to the destination satisfy 
frd > fSid, where i G {1, 2,..., M}. Second, we compare with 
the DBC scheme while setting g(j= 0, i.e., the destination can 
decode only if there is one node is transmitting per time slot. 
We refer to this scheme by the collision model of the DBC 
scheme (CM-DBC). This scheme is important to illustrate how 
the MPR capability at the destination is essential in the absence 
of sensing at the relay. Third, we compare with the SBC 
and DBC schemes but with /3i2=/32i=0 and 0 : 12 = 0 : 21 = 0 , 
respectively, to illustrate the effect of these probabilities and 
demonstrate the cases where their effect on the performance is 
essential. By setting /3i2=/?2i=0 and cti2=a2i=0, the relay 
can not transmit a relayed packet from Q ri by interfering on 
Sj , where i ^ j. 

We consider three cases for the channel conditions. 
In the first case, the system parameters are chosen as 
follows: P=W,R=l, p 2 Sud = 0.02, ^ 2>d =0.84, p^ 1)r =0.97, 
p 2 S2 r =0.93, and p 2 d =0.03. This case corresponds to an asym¬ 
metric channel situation, where the r-d channel and the si-d 
channel are weak while S 2 -d channel is strong. In Fig. |2] we 
plot the stability region for the cooperative schemes via vary¬ 
ing the value of w\ from zero to one by step 0.1. It is obvious 
from the figure that the SBC scheme significantly outperforms 
the CCMA scheme. The rationale behind this enhancement is 
that in the CCMA scheme the relay is restricted to transmit its 
packets in the idle time slots only. On the other hand, the SBC 
scheme adds to the relay the capability to, simultaneously, 
transmit its packets with the source terminal while controlling 
the interference probabilistically. This capability expands the 
stability region. 

It is worth noting that the two proposed schemes, SBC 
and DBC, can achieve exactly the same maximum stable 
throughput for S 2 , while this is not the case for Si. In the 
DBC scheme, the relay does not sense the channel, and hence, 
the relay may interfere with the source terminals. If the relay 
transmits a packet without sensing and interferes with S 2 , the 
destination with high probability decodes the packet of S 2 
first by treating the relay’s signal as noise, since the S 2 -d 
channel is strong, then decodes the relayed signal. This is 
not the case for si because if the relay interferes with si 
in the presence of these weak channels, s\-d and r-d, the 
destination most probably fails to decode both signals. The 
huge performance gap between the DBC scheme and the CM- 
DBC scheme demonstrates the importance of MPR capability 
at the destination in the absence of sensing capability at the 
relay. Removing the MPR capability causes a catastrophic 
reduction in the performance. 













Fig. 2: Stable throughput region for asymmetric channels 
configuration with weak relay-destination channel. 


Fig. 4: Stable throughput region for asymmetric channels 
configuration with strong relay-destination channel. 



Fig. 3: The maximum aggregate stable throughput for two 
source terminals in asymmetric channels configuration and 
weak relay-destination channel. 


In Fig. [3] we use the same system parameters as in the 
former figure. We plot the maximum overall stable throughput 
versus w±, which is the fraction of time allocated for Si. This 
figure shows that as wi increases the maximum aggregate 
stable throughput of the network decreases. The reason behind 
this is that the maximum stable throughput for si is equal to 
0.035 while that for s 2 is 0.9. Consequently, as we allocate 
more time slots for s 2 , the aggregate stable throughput of 
the network increases. Moreover, as shown in Fig. [2] the two 
proposed schemes, SBC and DBC, can achieve the same max¬ 
imum stable throughput for S 2 which dominates the aggregate 
throughput for any w\, and this is the reason why the aggregate 
throughput for the proposed schemes is close in Fig. [3] 

In the second case. Fig. [4] we plot the stability region 
for different cooperative schemes for the following system 
parameters: P=10,f?=l, p 2 Si d =0.8, Pg 2 (J =0.08, p 2 Sl r =0.85, 
p 2 2 r =0.9, and p 2 d =0.97. This case corresponds to an asym¬ 
metric channel situation, where the channels s±-d and r-d are 
strong while the channel S 2 -d is weak. It is important to note 
that the r-d channel is stronger than that in the former case. 


In this case both proposed cooperative schemes achieve the 
same stability region, and hence, sensing does not increase the 
stability region of the system. In the DBC scheme, the relay 
transmits according a random experiment. If the relay decides, 
randomly, to interfere with si, the destination can decode 
both signals because both channels, s\-d and r-d, are strong. 
Alternatively, if the relay interferes with S 2 , the destination 
most probably decodes the relay signal first, since r-d channel 
is strong, then decodes the source signal. Consequently, the 
MPR capability at the destination decreases the need to detect 
idle time slots, and hence, both proposed schemes can achieve 
the same stability region. Removing the MPR capability at the 
destination, in absence of sensing, yields to the performance 
of the CM-DBC scheme which is much worse than that of the 
DBC scheme. 

It is important to notice that both proposed schemes outper¬ 
form the CCMA scheme which restricts the relay to send only 
in the idle time slots. Another insight from this figure is the 
importance of pi 2 and P 21 in the SBC scheme and a: 12 and 
CK 21 in the DBC scheme to achieve this stability region. These 
probabilities have crucial effect in the asymmetric channel case 
as they enable the relay to utilize the time slots of one source to 
transmit the traffic of the other. Note that most of the packets 
of Si are transmitted directly to the destination due to the 
high gain of its direct channel. In contrast, S 2 suffers from 
low direct channel gain, so most of its packets are relayed. In 
the proposed schemes, the relay has the capability to interfere 
with Si, which has high direct channel gain, and send a relayed 
packet of S 2 , who suffers from low direct channel gain. This 
capability expands the stability region for both schemes. 

In Fig. [3 we plot the aggregate stable throughput of the two 
source network versus w\ . We use the same system parameters 
as in Fig. [4] It is obvious that the proposed schemes achieve 
higher aggregate throughput than that of the CCMA scheme. 
Unlike Fig. Q] the aggregate stable throughput increases as we 
allocate more time slots for si. In this case, the maximum 
stable throughput of si, Ag=0.98, is greater than that of 
S 2 , Ag=0.65. Consequently, as w\ increases the aggregate 
throughput of the network increases. 
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Fig. 5: The maximum aggregate stable throughput for two Fig. 7: The maximum aggregate stable throughput versus the 
source terminals in asymmetric channels configuration and transmission rate for symmetric channels configuration, 

strong relay-destination channel. 



Fig. 6: Stable throughput region for nearly symmetric 
channels configuration. 


In the third case. Fig. [6] we plot the stability regiorQ 
for a roughly symmetric channel configuration where both 
direct links, s\-d and s^-d, are strong. The system parameters 
are chosen as follows: p% d =0.75, p 2 2 d =0.8, p 2 Sl r =0.63, 
p 2 s2 r =0.73, and p 2 d = 0.85. The SBC scheme still exceeds the 
CCMA scheme, however, in this case the role of of /3i2 and 
/? 2 i diminishes due to the symmetric configuration. In this 
case, both source terminals have the same channel condition 
so the relay does not need to assist one source terminal at the 
expense of the other. Another insight from this figure is that, 
when the direct link channels are strong, we can achieve the 
performance of the CCMA scheme without even cooperation. 
However, the proposed schemes expand the stability region 
even when both direct links are strong. 

In Fig. [7] we illustrate how the proposed schemes respond 
when the channel capability to tolerate interference changes. 
We depict the effect of the transmission rate R on the maxi- 

2 We do not add the DBC scheme in the figure, because it is easy to realize 
from Fig. |d] that the performance of the DBC scheme in this case is exactly 
the same as that of the SBC scheme. 


mum aggregate stable throughput of the network for symmetric 
channel configuration. The parameters are chosen as follows: 
w=0.5, p 2 ud = 0.8, p 2 S2yd = 0.8, p 2 ur = 0.95, p 2 S2 r = 0.95, and 
p 2 (/ =0.96. It is clear that the stable throughput of the SBC and 
DBC schemes decrease slower than the CCMA and CM-DBC 
schemes. At low transmission rates, the channels can tolerate 
the interference, consequently, detecting the idle time slots is 
not important. It is obvious from the figure that the SBC and 
DBC schemes have the same performance for low transmission 
rates. As the transmission rate increases, the capability to 
sustain interference for all wireless channels decreases, and 
hence, it becomes essential for the relay to transmit only in the 
idle time slots. In the figure, as the transmission rate increases 
the performance of the two proposed schemes approaches to 
that of the CCMA scheme because the values of {Pij} 2 j—i 
and {ctij} 2 j—i begin to decrease to limit the negative effect 
of interference. As we increase the transmission rate more, 
the channels can not sustain any interference. The values of 
{Bij} 2 j_i ' n ^e SBC scheme diminish to be almost zeros, 
and the relay is restricted to send only in the idle time slots. 
This means that the SBC scheme boils down to the CCMA 
scheme. On the other hand, for the DBC scheme, as the trans¬ 
mission rate increases, the ability that the destination decodes 
the transmission of two nodes simultaneously decreases, and 
hence, the DBC scheme turns to the CM-DBC scheme. 

In Figs. 1 and E we illustrate the delay performance of 
the proposed schemes. First, in Fig. [8] we plot the minimum 
average delay encountered by the packets of S' 2 , under the 
constraint that A*=0.29. Moreover, we use the same system 
parameters as in Fig. [2] where the two proposed schemes do 
not achieve the same stability region. The maximum stable 
throughputs for s 2 , when A* =0.29, in the SBC and DBC 
schemes are 0.83 and 0.77, respectively. We do not plot 
the CCMA scheme to have a clear comparison between the 
plotted schemes, since the CCMA scheme performance is way 
worse than both. It is clear from the figure that the delay 
performance of the two proposed schemes in this case is close 
to each other for low A 2 . However, the SBC scheme delay 
performance slightly exceeds that of the DBC scheme. As A 2 
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Fig. 8: The average delay encountered by the packets of S 2 
for asymmetric channels configuration with weak 
relay-destination channel. 


Fig. 9: The average delay encountered by the packets of si 
for asymmetric channels configuration with strong 
relay-destination channel. 


increases the SBC scheme begins to significantly outperform 
the DBC scheme. It is obvious that the results obtained through 
queue simulation are very close to the results of the closed- 
form expressions derived in (fl9l > and ( 1301 ). This validates the 
soundness of the mathematical model. Moreover, the trade-off 
between the stable throughput and the average delay is clear 
where as the throughput increases the delay increases. 

In Fig. [9] we plot the minimum delay for si when Ag=0.81. 
We use the same system parameters as those in Fig. |4] where 
the two proposed schemes achieve the same stability region. 
In this case, the maximum stable throughput for S 2 in the 
SBC, DBC and CCMA schemes are 0.13, 0.13, and 0.038, 
respectively. The figure depicts that the proposed schemes sig¬ 
nificantly outperform the CCMA scheme. The rationale behind 
this is that the CCMA scheme allocates a large fraction of the 
time slots for S 2 to satisfy the constraint Ag=0.81, besides, 
the relay transmits only in the idle time slots. Consequently, 
the delay performance of the CCMA is much worse than that 
of the two proposed schemes. It is worth noting that even 
when the two proposed schemes can achieve the same stability 
region, the delay performance of the SBC scheme outperforms 
that of the DBC scheme. This is because the SBC scheme 
exploits the idle time slots to transmit the relayed packets, 
and the capability of the channel to sustain the interference to 
send simultaneously with the source terminal. Consequently, 
the SBC scheme exceeds the DBC scheme that exploits only 
the capability of the channels to tolerate interference. 

VI. Sensing at the relay versus MPR at the 

DESTINATION 

In this section, we illustrate how the MPR capability at 
the destination can compensate for the lack of sensing at the 
relay. Moreover, we derive the condition under which the two 
proposed schemes achieve exactly the same maximum stable 
throughput. A*. In this part, we assume that w\= 1, and hence, 
all time slots are allocated to si, i.e., we have only one source 
terminal si. 


Theorem 1. The maximum stable throughput for the two 
proposed schemes is given by 


a r = a-; 


1 Sl )(fs 1 d,+Ti)+ 1 Sl g r Sld 


1 ' Ti+g rd 

if the following condition is satisfied 

' 9 r Sld (T 1+ f rd ) 


frd - gld < min 


Tl + fsid 


(31) 


(32) 


~(Ti + 2 g s \) 


Ti(/, ld + Ti-^ id ) 

where 7i=(l- f Sld )f Bir . 

Proof of Theorem 1: In the SBC scheme, we can rewrite 
the stability condition in (ITOl) for «Ji=l as 


where 


A^min^r^f} 


p^ = (l-M(fs 1 d + T 1 ) + p 11 g: id 


u SBC = 

H'm 


frd 


-n SBC 
i Pi 


(1 — Pll)Tl + frd ~~ Pligrd 

The optimization problem in (IT2T > can be written as 

max mm{pl BC , p™ 0 } 

p ii 

subject to 0 < Pn < 1 

Lemma 1. If the following condition is satisfied 

f si ^ g r sA T i f rd ) 

frd - gfa < 


(33) 

(34) 

(35) 

(36) 

(37) 


Ti + f Sl d 

the solution of the optimization problem in ( D6D is given by 


As 1 * =(1-7 


Ti+9% 

Proof: See Appendix B 


Tl )(/, 1 d+Ti)+7-^ 5r <7s r i d (38) 


Tl+9rd 
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For the DBC scheme, we can rewrite the the optimization 
problem in ( f27l ) as follows 


min^rS 

subject to 0 < an < 1 


max 

an 


where 


Pi — (1 — aii)(/sid + T^+ang^d 
Q(ll frd 


u DBC = 




.DBC 


(1 - an)Ti + anfrd ~ «u g s rd 
Lemma 2. If the following condition is satisfied 

m + Q 2 (/ sld + ri) 


(39) 

(40) 

(41) 


frd gfct < 


Ti(f Sl d + Ti - g r sid ) 


-(T 1 + 2g s r f) (42) 


the maximum stable throughput for S! in the DBC scheme, 
which is the solution of the problem in ( 1391 ). is given by 


Af = (l- 7 


Tl+Prd 

Proof: See Appendix C 


1 Sl )(fsid+Ti)+ 1 s 1 9l 1 d 


Ti+g[ 


(43) 


rd 


The conditions ( l37l ) and ( 1421 ) establish the result in ( 1311 and 

<EDi. 


Theorem |T| states that if the values of f r d and gf d are 
close, i.e., they satisfy the condition in (f32l) . the two proposed 
schemes achieve exactly the same maximum throughput. The 
values f r d and g s r ' d depend on the variance of two channels; 
Si-d and r-d. From the definitions in ([2]) and (l44l) . the value 
of gf^ is close to that of f r d in two cases. First, if the two 
channels, s\-d and r-d, are strong, i.e., the destination can 
decode both signals with high probability. Second, if at least 
one of the channels, r-d or s\-d, is strong, i.e., the destination 
can decode the strong signal first then the weak one. 

We can easily map the insights from the obtained condition 
in ( l32l > to that in Fig. [2] and Fig. Q] First, in Fig. [3 the 
condition in (133 is violated for si, where we have two weak 
channels, s\-d and r-d. Hence, the maximum stable throughput 
of si in the SBC scheme is greater than that in the DBC 
scheme. Alternatively, in the same figure, the condition is 
satisfied for S 2 , where we have one strong channel S 2 -d. Thus, 
the maximum stable throughput of S 2 in the SBC scheme 
is exactly the same as that in the DBC scheme. Second, in 
Fig. m the condition is satisfied for both users, hence, the 
maximum stable throughput for si and S 2 is the same for the 
two proposed schemes. Moreover, the two schemes achieve 
the same stability region. From these results, we can realize 
that the MPR capability at the destination can compensate for 
the need for sensing to detect the idle time slots, if there is at 
least one strong channel to the destination. The strong channel 
facilitates the decoding at the destination, and mitigates the 
need of the relay to send only in empty channels. 

In Fig. [TO] we show numerically different channel variances 
that satisfy the condition in (132b . We plot the maximum stable 
throughput of Si versus f Sl d by varying <J 2 ld from zero to 
one for three fixed values a 2 d . The system parameters of 
this figure are chosen as follows: w\=\, P= 10, R= 1, and 
p 2 r =0.8. In the first case, where cr 2 d = 0.05, the two proposed 



Fig. 10: The effect of the channels, s\-d and r-d , on the 
maximum stable throughput of the two proposed schemes for 
u'i =1 for three fixed values for r-d channel variance. 


schemes achieve the same stable throughput approximately at 
/sid= 0 . 15 . As we increase cr 2 d , the performance of the two 
proposed schemes becomes closer than that of the former case. 
Ultimately, for the case of strong r-d channel of d =0.8, the two 
proposed schemes almost achieve the same performance for 
any cr 2 Sld , and this result emphasises the insights obtained from 
the condition in (133 . 

VII. Conclusion 

In this paper, we have proposed two cooperation schemes, 
and studied their impact at the medium access layer metrics 
such as stable throughput and average delay. In the SBC 
scheme, the relay senses the channel at the beginning of each 
time slot, and it decides either to transmit or receive packets 
depending on the sensing outcome. On the other hand, in the 
DBC scheme, the relay does not sense the channel, and it 
decides its operation in a random fashion. For each scheme, we 
derived the stability conditions for each queue in the system 
and characterized the stability region. Furthermore, we derived 
approximate expression for the average delay encountered by 
the packets. We illustrated how the SBC scheme significantly 
outperforms existing cooperative schemes. The SBC scheme 
exploits the available resources more efficiently than other 
cooperative schemes because the relay not only utilizes the 
idle time slots, but also interferes with the source terminals in 
a mild way to mitigate the adverse effects of interference. 

Moreover, we demonstrated that the MPR capability at the 
destination can compensate for the relay need to sense the 
channel. Although the relay in the DBC scheme does not 
sense the channel, our results show that the DBC scheme can 
achieve, under a certain condition, the same stability region as 
that of the SBC scheme. 

Appendix A 
Derivation of g 

The term denotes the probability that the link ( m,n) 
is not in outage in presence of interference from node I. Thus, 
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we can express it as follows 


9mn Vmn 4“ (^ Lrm ) ^ rnn (44) 

where v I mn is the probability that the node n £ L successfully 
decodes both packets transmitted from in £ T and I £ T, 
on the other hand, h I mn is the probability that the node n 
successfully decodes the packet transmitted from m by treating 
I as noise. Let X and Y be two independent exponential 
random variables with parameters 71 and 72 , respectively, and 
p[x) and p(y) be their probability density functions. We define 
two deterministic variables i] and r/j. To derive the expression 
of we first define the region 


Vi) = {(£, y) : X > 17 n y > (lx+y >771} ( 45 ) 

= {(x, y) : x > 77 n y > max[rj , ??i - x]} 

= {(x, y) : 771 - 77 > x > 77 (~l y > 771 - x} 

U {(x, y) : x > 771 - 77 n y > 77}, 

then, we figure the following integration 


^( 77 , 771 , 71 , 72 )= // p(x, y)dxdy 


$ft(77,77i) 

r/i—ri 00 


p(x) / p(y) dydx 


T) Tj — X 

OO OO 


p(x) / p(y) dydx 


(46) 


Vl-V V 

o72J?l 


7i e 


-(exp(-(7i-7 2 )r;)-exp(—( 7 !—7 2 )(tji— 77 ))) 


71-72 
+ exp(77i7i + 77 ( 72 - 71 )) 


Hence, we have 


= V{ 


2 R -1 2 2R —1 1 1 


P ’ P ’ p 2 ’ " 2 


’ p 2 

m,n r/,n 


(47) 


To derive the expression of /i^ n , we first perform the 
following integration 


^(77,71,72) = P{i? < log(l + 


Px 


= / exp 1 


Pt/ + 1 
(2*-1)(P7/ + 1) 7 i 


P 


)} (48) 

)p{y)dy 


72 exp (-7771) 


72 + P7177 

Hence, the probability h I mn is given by 


— 1 1 1 

^mn 7) , —n , —9 ) 


Pm,n Pl,n 


(49) 


Appendix B 
Proof of Lemma 1 

Taking the derivative of /if BC and with respect to /3u 
yields the following 

dpf BC 


d/3n 


= 9 r s ld ~fs 1 d~T 1 


dp S * C / rd K d -Ti-/ Sld )(Ti+/ rd -^ii(Ti+^)) 

00 n 


(50) 

( 51 ) 

((l-0n)Ti+ /„,-&!<$) 

| fUTi + g^)m+/ Sld +/3n(<4 d + Ti ~ fs ld )) 

((l-/3n)Ti + / rd -/3n^) 2 

From (l50t . we can see that /r) BC is a monotonic ally decreasing 
function in (3u because, from definition, f Sld is greater than 
Slid- 11 11 4ic other side, from ( ITU , we can show that /i5 BC is 
monotonically increasing in (3u if the following condition is 
satisfied 

f Si ^ 3sid(-Pl + frd) 

Jrd 9 r d — 


7"l + fsid 


(52) 


Since the value of p,f BC at /?n=0 is greater than the value 
of /r* BC at /?n =0 and /if BC decreases monotonically with 
/3n while /r® BC increases when (f52l) is satisfied, the optimal 
solution of (l36l) occurs when /r'j iBC = /r^ BC . Thus, under the 
above condition in (l52l) . the optimum 3\-\ is given by 


Pn = 


T\ 


(53) 


Ti + a? d 

and, hence, the maximum stable throughput for si is given by 


As 1 * = ( 1 - 


Ti 


T l+9 rd 


-)(fs ld +Tl) + 


Ti 


Ti+ 9 s r y 


(54) 


Appendix C 
Proof of Lemma 2 

It is clear, from ®, that Pi BC decreases monotonically 
with an. On the other hand, the derivative of ti° bc with 
respect to an is given by 

dp™ C _ MT 1 + f 3ld +2 ail (g^ d -T 1 -f aid ))(T 1 + ail (frd-T 1 -g a r 1 d )) 
<9an ((1 —aii)Ti + an/ rti —/3ng“J) 2 

frd(frd—Ti—g^)(aii(f sl d+T 1 )+a'li(g^ l d — fs 1 d. — Ti)) 

((l-QiilTi+aii/r-d-ang ^) 2 

(55) 

It is obvious from (TtII) that p^ BC is a quadratic over linear 
function in an, and hence, it is a concave function in an 
for the range of an from zero to one l20l . Besides, /i|) BC 
increases at an= 0 , because the slope of p,° BC at an =0 is 
positive and given by 

_ frdTljTl + /sid)_ 

((1 - au)Ti + an/rd - anSrd) 2 

Note that the two functions ti![ )BC and /r° BC satisfy the fol¬ 
lowing three conditions. First, the value of // l ) Bf: at an = 0 
is greater than the value of p,^ BC at the same point. Second, 
// l | )B<: is monotonically decreasing function in an. Third, p% BC 
is concave function with positive slope at an=0. Therefore, 
the optimum value of an is obtained at the intersection point 
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between and /z™ c if the slope of is positive at this 
point because this means that /t™ c increases monotonically 
with an from an= 0 until the intersection point. has 

a positive slope at the intersection point if the following 
condition satisfied 


frd g r d < 


m+g^) 2 (/ si d+ri) 

Ti(f Sl d + T\— g r Sl(] ) 


~(.T 1 + 2g° 1 d ) (57) 


Under the above condition, the optimum solution of 
occurs when // l l )BC =/i,'’| iC , hence the optimum value of an 
is given by 

' (58) 


ON i — 


Ti + g'A 

and, the maximum stable throughput for si is given by 

, 1 Sl )(fs 1 d+Ti)+ 1 s i gl 1 d 


Ti+g; 


rd 


Ti+g[ 


(59) 


rd 


which is exactly the same stable throughput for the SBC 
scheme in (l54l >. 
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